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A NOVEL, AXISYMMETRIC, ELECTROSTATIC COLLECTOR 
FOR LINEAR BEAM MICROWAVE TUBES 
by Henry G. Kosmahl 
Lewis Research Center 

SUMMARY 

An analytical model of a novel, axisymmetric, electrostatic depressed collector for 
linear beam microwave tubes is introduced. The model features schemes for independ- 
ently effecting the sorting of low- and high-energy electron classes, prevention of back- 
streaming secondaries and primaries, and a geometry that makes the size of the beam 
unimportant for the shape of trajectories. 

The boundary -value problem is solved by developing a Green's function in elliptical 
harmonics. The trajectories are derived from Lagrangian equations including accurate 
space-charge effects. The magnetic boundary -value problem is solved for an ideal 
shield with an aperture for the exit of the spent beam. The charge conservation is 
treated for two assumed velocity distributions across the beam: a nearly parallel multi- 
velocity beam of constant current density and a diverging multivelocity beam . Examples 
of digitally computed electron trajectories are shown. The performance of the collector 
and the problem of optimizing the overall tube -collector efficiency is discussed. 


INTRODUCTION 

The residual energy in spent electron beams emerging from the exits of microwave 
tubes represents not only a loss in efficiency but also an unnecessary heat dissipation 
burden. The advent of space communication systems with their demands for highest 
realizable overall efficiencies and, if possible, for the elimination of the cooling problem 
makes the recovery of this energy mandatory. Since we deal here with kinetic energies 
of electrons, an obvious approach is to slow down the electrons and collect them on sur- 
faces at the lowest possible potentials. The name "depressed collector" appropriately 
describes devices fulfillir^ this purpose. 

Depressed collectors have been subject to numerous experimental investigations and 


a few simplified and quite superficial analytical treatments since the late fifties. 
Wolkstein (ref. 1) and Sterzer (ref. 2) pioneered in successfully applying a one- or two- 
segment depressed collector, respectively, to low -efficiency traveling wave tubes. 

Their collectors consisted essentially of one or two segments of cylindrical tubings 
insulated from each other and shielded from the focusing magnetic field. Since in low- 
power, low -efficiency tubes both the radial and axial velocity spreads are small, a 
significant improvement factor in efficiency in comparison to "undepressed” operations 
was achieved. The degree of depression, while being far less than a theoretical maxi- 
mum, was determined by the onset of backstreaming electrons, either reflected pri- 
maries or secondaries. Moreover, since the metallic boundaries of the segments did 
not coincide with natural equipotential surfaces, strong fringing fields ("lens effects”) 
were formed. It is obvious that the presence of such fringing, strongly curving fields 
will, in general, prevent the collection of electrons at lowest possible potentials and will 
also cause the backstreaming of many electrons, not necessarily due to a lack of energy 
but rather because of unfavorable deflections. 

Attempts to improve the overall tube efficiency by collector depression in high- 
power tubes with large velocity spreads were many, but not much success was achieved 
with conventional techniques. With a few exceptions, the results were published only in 
internal reports. 

New approaches which eliminate some of the deficiencies of the earlier works were 
introduced by Sauseng (refs. 3 and 4) and Branch and Mihran (ref. 5) and Neugebauer. An 
effective method to suppress secondaries and their backstreaming was introduced at 
Bell Telephone Laboratories, and later at Hughes by Sauseng, by collecting electrons in 
negative electric fields. To enhance the velocity sorting, Sauseng (ref. 4) employed an 
asymmetric magnetic deflection. This method, although it is very effective in sorting 
off-axis as well as on-axis electrons, seems to have a basic deficiency due to unequal 
action of the magnetic field on electrons parallel and perpendicular to it. 

Branch and Mihran (ref. 5) and Neugebauer found the interesting property of elec- 
trons deflected by a long symmetric negative spike to "focus" at one point after coming 
down from the apex of the trajectory. 

The approach pursued in the present investigation takes a basic configuration of a 
sphere with a central cone and a short protruding spike as a flexible geometry to perform 
three functions effectively at the same time; that is, (1) to velocity-sort the electrons, 

(2) to slow them down to minimum realizable velocity at collection, and (3) to prevent 
the backstreaming of secondaries and reflected primaries into the interaction region of 
the tube. 
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METHODS OF SOLUTION 


The solution of the potential problem is derived in appendix A with the powerful and 
elegant tool of Green's function. 

As shown in figure 1, the cone, the spike, and a part of the spherical surface are 
assumed to be at the potential -1.5 V^, the surface adjacent to the injection hole at V = 0 
and a continuous distribution from V = 0 to -1. 5 over the rest of the surface. (All 
symbols are defined in appendix F. ) Here, the cathode is at a potential -1. 0 V^, the 
tube output at V = 0. In the tube considered for future applications, only very few elec- 
trons will have energies larger than 1. 5 V^. Thus, there will be little interception 
on the spike and cone. The distribution of the potential between hole and cone on the 
boundary of the sphere is a variable to be chosen such that sufficient sorting is obtained 
with the highest possible energy recovery. It is obvious that a purely electrostatic de- 
flection requires some small part of axial energy to be converted into radial components 
which produce the sorting. 

It is apparent from reviewing previous efforts with "conventional" collectors that 
their failure to perform effectively the first and third required functions limited their 
usefulness in applications to more efficient tubes. In contrast to this deficiency, an 
inspection of the equipotential lines and collecting electrodes, as well as of the apertures 
for the electron passage, shows the effectiveness of the present design to fulfill all re- 
quirements. Since the equipotential lines slope upward, the direction of force on elec- 
trons is away from the axis of symmetry. Also, the sharp spike produces an indentation 
in the equipotential lines in the immediate vicinity of the axis. Thus, slight perturba- 
tions such as noise, small plasma fluctuations, or residual signals will cause even elec- 
trons going straight up to be quickly deflected away from the axis, where they become 
subject to sorting and collection. These processes are discussed in the section Analyti- 
cal Results. 


Initial Beam Conditions 

Any exact calculation of electron trajectories in the collector region, which is 
treated in appendix B, presupposes a knowledge of the electron vector velocity distribu- 
tion in the entrance plane as a function of time, that is, for all phases of the full radio- 
frequency cycle. If azimuthal symmetry is assumed, the problem is reduced to the 
knowledge of a distribution function f^(*, z,*", z, (p,t) in a six-dimensional phase space, 
the subscript i designating the initial injection plane. Since presently available non- 
linear (large signal) computer programs are one -dimensional, fj is not yet known. How- 
ever, from fairly accurate nonlinear three-dimensional calculations (unpublished current 
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work of this author), it is evident that most electrons enter the collector at angles a. 
between 0° and 15° at radiofrequency voltage levels from zero to the degree of de- 
flections being determined almost entirely by the distance from the axis rather than by 
phase or radiofrequency voltage. Since these calculations were carried out for, among 
others, a confined flow through a magnetic field of magnitude two to three times the 
Brillouin field Bg, terminated by an ideal shield, it is obvious that the magnetic deflec- 
tion largely dominates the problem. Thus, in absence of exact data, a deflection angle 
proportional to the radial position of the electron may be assumed. Aside from the 
direction, the magnitude of the velocity is of interest. Especially at high frequencies 
(e.g. , larger than X-band), the beam diameter is so small compared to collector dimen- 
sions that the electrons may be considered as coming from a point source, but with an 
equivalent distribution ranging from about zero to 2 V^, depending on radiofrequency 
phase. Another important feature of the beam is its degree of laminarity. The problem 
is complicated by the number of possible cases, including klystrons and traveling wave 
tubes (TWT) and solenoid and ppm focusing. In a well -focused, confined-flow klystron 
beam through a solenoid, good laminarity may be assumed everywhere, except in the 
output cavity during decelerating fields. In a solenoid -focused, coupled -cavity TWT, the 
onset of nonlaminarities may begin two or three cavities prior to the output. In ppm- 
focused tubes, nonlaminar behavior may prevail over the end of the tube due to strong 
scalloping. Thus, an exact analysis of ppm -focused traveling wave tubes is so complex 
that efforts toward providing it should be critically evaluated with regard to its impor- 
tance. 

Fortunately, all indications show that the majority of the electrons will be contained 
within a small angle in all cases except electrostatically focused tubes. Thus, their 
radial velocities can be reduced, though not eliminated, sufficiently to permit a substan- 
tial energy recovery in the collector. Another argument to be accepted is the fact that 
even if the exact electron distribution were known, no single or simple enough refocusing 
method would eliminate the radial spread for all electrons at all phases. 


Space-Charge Effects 

It is appropriate at this time to evaluate qualitatively the effects of space charge on 
the operation and performance of the collector. Whereas at low frequencies (f < 2 GHz) 
the current density at the entrance into the collector is of the order of 1/2 to 3 amperes 
per square centimeter, this number increases to approximately 40 to 100 amperes per 
square centimeter at 12 gigahertz, depending on perveance. An examination of the 
trajectories and equipotential lines obtained for the two frequencies and presented in 
figures 1, 2, 4, and 3, respectively, demonstrates clearly the highly undesirable effects 
of high space charge on the trajectories and performance. In the low-frequency design 
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(a) Electron injection radius, 0.4 centimeter. (b) Electron injection radius, 0. 8 centimeter. 

Figure 2. - Trajectories and equipotentials for a collector with a 0. 1 R long spike. Magnetic field intensity, 0; current 
density, 0.5 ampere per square centimeter; injection angles. O'". 5°, W. 
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(a) Current density, 30 amperes per square centimeter. 
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potential 
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180" 

(b) Current density, 90 amperes per square centimeter. 


Figure 3. - Trajectories and equ (potentials for a collector without a spike. Magnetic field intensity, 0; electron injection 
angles between T and 16"; injection radius. 0. 1 centimeter. 
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Electrode 



180° 180° 

(a) i\onlnverted cone. (b) Inverted cone. 


Figure 4. -- Trajectories and equipotent'als for a magnetically confined flow beam and Brillouin flow. Frequency, 
2 qig?.hert7; current density, 0.5 ampere per square centimeter. 


(small space charge) shown in figures 1, 2, and 4, the Poisson solution hardly differs 
from the Laplace solution; that is, the external fields dominate the potential shape and 
electron flow. On the other hand, in the high-frequency case (fig. 3), not only are the 
space -charge fields much stronger (in the vicinity of the injection hole) than the space - 
charge -free fields produced by external potentials, but there is an acute danger of space 
charge blocking the flow of electrons and a buildup of a virtual cathode. The deformation 
and depression of equipotential lines causes an excessive deflection of electrons radially 
outward, thus decreasir^ the degree of energy recovery . In addition, the direction of 
the flow is clearly space -charge dominated, making a planned velocity -sorting very dif- 
ficult and very sensitive to space -charge fluctuations and possible instabilities. In addi- 
tion, a refocusing of such high-density beams is not practical. We must, therefore, 
conclude that a substantial dilution of high current densities down to 1 -ampere -per - 
square -centimeter levels is required if collectors designed for h^h -frequency tubes are 
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to be made feasible and approximately as efficient as low-frequency designs, A possible 
arrangement to accomplish this task is shown in figure 5 and is discussed in the section 
Analytical Results. For the case of the highly desirable small space charge in the col- 
lector region, even large errors in computing will not change the trajectories sig- 
nificantly. This, in turn, permits the application of simplified methods for obtaining an 
expression for p^. For this reason, a semi-quantitative approximation is made regard- 
ing the velocity vector distribution across the beam. This approach is treated in 
appendix C in connection with the continuity equation. 

The leakage of the magnetic focusing field into the collector region through the 
aperture in the shield is solved in appendix D. An existing axisymmetric expansion of 
the vector potential is used within the diameter of the solenoid, using the field along the 
axis. In addition, a solution valid everywhere is derived by transforming the scalar 
potential, resulting from solving the magnetic -boundary -value problem on the shield and 
in the aperture from spheroidal oblate coordinates into spherical coordinates. 

The relations relating tube efficiency, collector efficiency, and overall efficiency 
are discussed in the section Analytical Results. 

The international system of units is used throughout. 
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DISCUSSION OF COLLECTOR OPERATION 


General Considerations 

Consider the computed electron trajectories in figures 1 to 4. They represent the 
motion of a spent electron beam which, while being axisymmetric, has a wide distribu- 
tion in energy, radial and axial velocities, and radial distance from the axis at the time 
of entry into the collector. It is clear, therefore, that the task to collect all electrons 
which started within a given range of initial kinetic energy at close to zero final energy 
on the same electrode regardless of angle, radial position, and radial velocity is an 
impossible proposition for a static arrangement of electrodes. This difficulty can, how- 
ever, be somewhat mitigated by making the beam radius at entry small compared to the 
distance the electrons have to travel from injection point to collection point. In this case 
the electrons appear to be coming from a point source, thus rendering the distribution in 
radius unimportant. For the range of remaining variables (i.e. , the distribution in 
energy and direction), a compromise must be made such that the impact kinetic energy 
is a minimum. 

Another important design feature is the degree of flexibility in accomplishing the 
sorting into energy (velocity) classes. Some axisymmetric collectors (ref. 5) use a 
long, protruding spike as the only sorting (deflecting) element. In the present configura- 
tion the sorting is effected separately by a relatively short spike, which acts mainly on 
the more energetic electrons, and by the potential distribution on the sphere and the cone, 
which determines mainly the slope of the equipotentials in the vicinity of the entrance. 

The advantage of this flexibility becomes apparent if particular velocity classes are in- 
jected in a preferential direction which requires a different degree of sorting (deflecting 
force) for different energy classes. For high collector efficiency, only as little radial 
deflection should be applied as is necessary for steering the electrons of a given energy 
class to an electrode, where they can be collected with a minimum of total kinetic energy, 
that is, axial and radial components. 

Still another important requirement concerns the prevention of backstreaming of 
primaries and secondaries. As mentioned earlier, streaming of secondaries may be 
eliminated entirely if they are generated on the "upper" side of the electrodes facing 
negative electric fields. The analysis and visual inspection of trajectories shows that a 
great majority of trajectories are collected on the upper side of the plates (except on the 
cone and the spike). However, a small fraction of electrons will strike the edge of the 
underside causing the generated secondaries to stream to the next higher potential elec- 
trode. In order to keep the area of the struck edge as small as possible, the electrodes 
should slope "upward" away from the axis. 

As far as the backstreaming of reflected primaries is concerned, the best way of 
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keeping this effect to a minimum is by designing the diameter of the apertures in the 
electrodes as small as possible. The limit on the size is dictated by the onset of pri- 
maries impinging on the undersides, which prevents energetic electrons from being col- 
lected at yet lower potentials. Theoretically speaking, there is, in the axisymmetric 
configuration discussed herein, no possibility of preventing backstreaming of primaries 
and secondaries which move down exactly along the axis since there is no radial (de- 
flecting) force along it. However, the number of electrons within a given radius is pro- 
2 

portional to r , where r is the distance from the axis. Thus, the percentage of elec- 
trons traveling close to and parallel to the axis is small. In addition, space-charge 
effects, radiofrequency instabilities, and random motions will tend to deflect electrons 
from the axis. Nevertheless, the axis of symmetry is the "weakest” spot and the only 
known deficiency of the electrostatic, axisymmetric collector and should be given proper 
attention. 


Analytical Results 

Figures 1 to 3 show electron trajectories for various beam parameters and collector 
configurations in the absence of magnetic focusing fields (i.e. , B = 0). Figure 4(a) shows 
trajectories for magnetically focused beams which enter the collector after passing 
through a hole in an ideal magnetic shield p = °°). All confined-flow focused 

beams experience a deflection during the transition from B = into B = 0, which is 
due to the flux at the cathode. This "magnetic" deflection is absent for Brillouin- 
focused (ii/„ = 0) and electrostatically focused beams. Although the rotational (magnetic) 

c 2 2 *2 

energy (m/2)* (sin ^’))cp at the apex is always small compared to the initial injection 
energy (m/2)* a is large (e. g. , ::?l/3 e^) and the collector efficiency, particularly 
for the low-velocity classes, is relatively small for the equipotential configuration shown 
in figure 4(a). This deficiency can be much reduced, as shown in figure 4(b), by making 
the equipotential s more perpendicular to the trajectories at the apex. 

It is important to shield the collector region (as perfectly as feasible) from the main 
focusing field for two reasons; First, the amount of rotational energy which cannot be 
recovered at impact is reduced. And second, "back focusing" of reflected primaries and 
generated secondaries into the interaction region of the tube along magnetic flux lines can 
be prevented. 

In figures 1, 2, and 4 the space -charge density is small and affects the trajectories, 
compared to space -charge -free cases, only very little. In figure 3, computed for injec- 
tion current densities of 30 and 90 amperes per square centimeter, the space-charge 
forces dominate the trajectories in the vicinity of the injection hole. Notice the heavy 
depression of equipotential s in the injection region. The spike has been removed en- 
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tirely in order not to obscure the space -charge effects. It may be seen that space charge 
provides too high a degree of sortii^ for low-energy classes because approximately one- 
half of the axial energy is transferred into radial energy, with a subsequent severe loss 
in collector efficiency. Notice further in figures 3(a) and (b) that the apex of the trajec- 
tory for energy class 0.2 is only slightly above the - 0. 1 equipotential line, yielding a 
collection efficiency for this class of approximately only 50 percent. In contrast to this 
case, the low -space -charge cases in figures 1, 2, and 4 show the energy at the apexes to 
approach 80 to 95 percent of their respective injection energy, thus permitting in princi- 
ple a high collector efficiency. In addition, it is very difficult to make a planned 
velocity-sorting efficient in cases of space -charge -dominated flow. There is a danger of 
virtual cathode buildup and subsequent space -charge blocking and possible breakout of 
oscillations due to instabilities and fluctuations in high-space -charge densities. We must, 
therefore, conclude that a substantial dilution of high current densities down to a level of 
a few amperes per square centimeter or less at the injection hole is required if collector 

designs for higher -frequency tubes are to become feasible. 

A possible arrangement to accomplish this task is shown in figure 5: the beam is 

permitted to expand in diameter by a factor of 5 or more, that is, the current density 
is decreased by a factor of 25 or more to the 1 -ampere -per -square -centimeter level. A 
simple electrostatic lens compensates the outward radial electron velocity +Uj. by pro- 
ducii^ an inward component -u^. on an electron located at radius r in the plane of the 
lens. Had all electrons exactly an outward radial velocity u^. which is strictly propor- 
tional to r, then an exact compensation of radial velocities would be possible and a 
close -to -parallel flow would enter the collector. Unfortunately, in a real spent electron 
stream, many electrons have radial velocities which are not strictly proportional to r 
and also the trajectories have positive and negative slopes near the lens due to a wide 
spread in velocities. Therefore, it v/ill be possible, at best, to reduce, on the average, 
only a fraction of the radial spread, with a resulting substantial improvement in collector 
efficiency compared to entirely uncompensated cases. A magnetic lens permitting an 
oscillation -free adiabatic expansion to a larger Brillouin radius could be used for this 
purpose, as well as a combination of electric and magnetic fields. However, a more de- 
tailed discussion requires a precise knowledge of three-dimensional trajectories. 


Relations Between Tube Efficiency, Collector Efficiency, 
and Overall Efficiency 

A number of efficiency definitions are possible and needed by the tube designer and 
system engineer using the tube. Consider f^re 6 on which the overall efficiency 77 ^^ 
is plotted against the coUector efficiency with tube efficiency Tjj. as parameter. In 
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Figure 6. - Relation between overall tube efficiency and 
collector efficiency. ' “^c*] * 


this plot, circuit losses and interception losses prior to the collector are neglected. K 
is the beam power I^V^, then 
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A more accurate formula is given in equation (4) 

% 



( 2 ) 


(3) 


(4) 


where 77 ^, and P.j^^ designate the circuit efficiency, electronic power conversion 

efficiency, and the power intercepted prior to entering the collector, respectively; 

Pint is diff icult to determine but it can be written as (e) /Iq, where <e)isthe 
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normalized average volt energy of the intercepted electrons and 1 ^^^^ is the intercepted 
current. At saturation (e) is usually between 0. 1 and 0. 8 , depending on particular 
circumstances. 

In equation (3) the summation is over n collector electrodes, each collecting a 
current Ij^ at a potential Vj^ below the potential zero of the output section of the tube. 
Thus, 2 Ijj ■ is the kinetic power in the spent beam recovered by the collector and 
(1 - 7 ?()Pq is the kinetic beam power entering the collector. Note that is less than 
the electronic efficiency of the tube tj^. 

How does one optimize the only efficiency of interest to the system designer ? 

Although a quantitative answer is not yet available, a quantitative estimate can be easily 

produced. The electronic tube efficiency 77 is proportional to the square of the total 

2 ^ 

beam coupling coefficient k times the time average of the product ij^(t)vj^(t) cos (p(t), 
where ij(t) and Vj(t) are the time -dependent amplitudes of the radiofrequency convec- 
tion current and radiofrequency voltage, respectively, and ^(t) is the phase between 
them. Thus, 


Vq ~ ii(t)v^(t) cos (p{t) 

The theoretical maximum of i^ is 2 1^. In well-designed klystrons i^ approaches 
values of 1. 6 to 1 . 8 1^, and in traveling wave tubes 1. 2 to 1. 4 1^. The maximum of the 
radiofrequency output voltage Vj^ must be smaller than that which causes the reflection 
of the slowest electron in the output gap (i. e. , Vj^ < 1. 1 Vq). Although the adjustment of 
the output voltage to such a high value may produce an optimum of 77 ^, it certainly will 
not yield the highest possible 77 ^^. It is clear that the spent beam will contain many elec- 
trons with energies close to zero whose residual kinetic energy can not be recovered and 
which only cause high interception at undepressed potentials. In addition, such a stream 
will have a very high degree of axial and radial velocity spread, which makes effective 
sorting less efficient. We have to conclude, therefore, that an optimum in klystrons 
will probably occur at output radiofrequency voltages slightly less than those producing a 
maximum in 77 ^. A quantitative discussion would require the knowledge of exact three- 
dimensional trajectories. 

In traveling wave tubes the situation is different. Because the TWT circuit imped- 
ance Z is smaller than in klystrons, the radiofrequency voltages v-, are almost always 

I 1 3 

well below (e.g. ^ Vq), and the slowest electrons in the stream have still a 

significant kinetic energy. It appears, therefore, that depressed collectors should have 
a better potential of improving TWT efficiency than that of klystrons. 
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Collector Losses 


It is appropriate at the conclusion of this discussion to review basic limitations and 
losses of this collector. Were exact velocity vectors known, it would then be possible 
to give an accurate account. In absence of such knowledge, only a useful estimate can 
be given. 

An electrostatic collector requires a certain energy, fortunately only a small 
amount, to accomplish velocity -sorting. The amount of this energy which represents a 
basic loss of efficiency can be measured directly as the difference between the volt 
energy at the apex and the initial volt energy at the start. For small angular spread 
this difference can be made small, perhaps 5 to 10 percent. 

Another source of losses stems from the finite number of collector plates. For a 
large number of plates the electrons will be losing roughly one-half of the potential dif- 
ference between two adjacent stages. For a five -stage collector, this amounts to approx- 
imately 10 percent loss in efficiency, for a 20-stage collector only 2. 5 percent loss. An 
additional loss results if there are electrons with energies larger than where 

Vcath collector plates are below cathode potential. This loss is very signifi- 

cant in klystrons and less important in TWT's and may range between 5 and 25 percent. 
Still other losses may occur due to secondaries and ions, if present. 

All the aforementioned deficiencies can be reduced to a tolerable minimum with 
careful design and effort. There is, however, still present a different source of defi- 
ciency which is much more difficult to control and eliminate than the aforementioned 
sources. This deficiency is related to inadequate velocity -sorting caused by large 
angular spreads and can be only partially eliminated with refocusing. In figure 1(b), 
for example, the electron of energy class 0.2 and starting with 6-^ v/ill easily land on 
plate V = -0. 15 Vq. The same electron in figure 1(c) starting at 3*^ will just miss the 
-0. 15 Vq plate. Similar considerations apply to other trajectories. It is estimated that 
velocity-sorting losses can be reduced to 5 percent with refocusing. 

It therefore appears feasible that a carefully designed collector with low space charge, 
a refocused beam, and having 10 or more stages may approach 80 to 90 percent recovery 
efficiency from a spent beam. This v/ould represent a drastic improvement over older 
designs. 

Figure 7 shows an artist's concept of the collector exposing the apertures and elec- 
tron trajectories. Figure 8 shows the actual collector mounted on top of an S-band 
klystron amplifier. The individual electrodes are mounted in a rack and held separated 
by insulators. 
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CONCLUSIONS 


A model of an axisymmetric, electrostatic depressed collector for linear -beam- 
type microwave tubes is given. The model eliminates to a high degree the deficiencies 
of ’’older" depressed collectors and shows promise to work effectively even when applied 
to a spent beam with a large radial and axial velocity spread. 

It is concluded that the geometry of the beam and the collector must be chosen such 
as to render space -charge effects relatively unimportant bn trajectories and collector 
efficiency. In confined-flow focused beams, all electrons experience a deflection which 
results in a substantial loss in collector recovery efficiency for low-energy electrons 
and a small loss for higher energies. It is concluded that as small a ratio of the focus- 
ing field to the Brillouin field B/Bg as feasible should be used to enhance the collector 
efficiency. Such losses do not occur in focusing schemes where the magnetic field at the 
cathode is zero or small or when the beam has been refocused. A qualitative argument 
is derived for achieving an overall optimum efficiency for the tube -collector combin- 
ation. A significant improvement in the state -of-art of tube -collector efficiency is be- 
lieved to result with the optimization of the designs, with collector efficiencies approach- 
ing 90 percent for a less -than- 10 -electrode configuration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 16, 1970, 

164-21. 
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APPENDIX A 


SOLUTION OF THE BOUNDARY -VALUE PROBLEM 
WITH GREEN'S FUNCTION 

Common electrical boundaries, such as a spark gap or a spike, can be well approxi- 
mated by an axisymmetric conductor of the shape of a prolate ellipsoid of rotation. In 
the case of a spike approaching the shape of a sharp needle, the smaller semi-axis b 
of the prolate spheroid (e. g. , an axisymmetric ellipsoid) may be shrunk to the limit zero 
or to any other fraction b/c < 1, where c is the larger semi-axis. 

We solve the boundary -value problem with the method of Green's function. After the 
Green's function appropriate to the selected geometry has been found, the determination 
of the potential is reduced to a quadrature by means of the equation 


V(P) = y* G(P,Pq)-P^(Pq) d-^ + f Vg(S)^^Jj dS (Al) 

volume "surface S 

Here G(P, P^) is the Green's function of reference points P and source points P^, 
is the space -charge density inside the volume bounded by the surface S along 
which the potential is Vg. 

To derive the Green's function for our basic configuration of a cone with an included 
half -angle 9, a spike at its apex represented by the two semi-axes c and b and sur- 
rounded by a sphere of radius R (see fig. 9), we use the coordinate system rj, and (p 
of spheroids of rotation. The coordinates ri and | are related with the coordinates z 
and r in the cylindrical coordinate system by the equations (ref. 6, ch. 5.28): 
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z 

t 



Figure 9. “ Derivation of Green's function in spheroidal (elliptical) harmonics. 


Substituting equations (A3) and (A4) into equation (A5) yields 


or 


cos ^ 

J 2 ^ ,2 T 
- 1 


^ = cos B 





(A6) 


Putting cos B = Constant = jj.^ into equation (A5) we obtain on the surface of the cone 
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The difference between and is largest for = 1 and diminishes to zero as 
r] becomes large. Thus deviations between and occur only at the joint of the 
spike and the cone, where no trajectories are present. A hyperboloid 


, 2.2 




= 1 


could be used exactly in place of the cone to preset the inclination of the equipotentials 
(see fig. 9 ), the only reservation being that a curved surface is more difficult to produce 
than a cone. To obtain the expression for the Green's function, write the general solu- 
tion of Laplace's equation in coordinates of a prolate spheroid; 


with 


^ ^ ^ ~mn l^mn^m 

(A7) 


(A 8a) 

«mn = A*P™(r7) . B’Q-(t 7) 

(A 8b) 

= c cos nup + D sin nup 

(A 8c) 


Since we seek a solution which is finite everywhere, including the origin and infinity, we 

must reject Q™(^) inside because lim Q(|) = «>, and V{rf} outside since lim P(r]) = 

^—1 77—00 

The range of the independent variables 77 and ^ appearing in the formulation of equa- 
tion ( 7 ) is 1 < 77 < 00; -1 < ^ + 1 . Because in the presence of the spike the line 77 = 1 

is excluded, Q(t 7) must be admitted. Since our problem involves conical boundaries, the 
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solution (eqs. (7) and (8)) must be expanded in terms of harmonics with noninteger val- 
ues n so chosen that the P^(|) are zero on the cones. That is, 
satisfying the equipotential requirement on the surface of the cone, cos t 9 = ju « | = 
Constant, Since for noninteger values n Pjj(l) has a logarithmic singularity at | = +1, 
we must use the function Pj^(-|), which is regular in the region 5 > 0, where Green’s 
function is derived. 

To solve this problem, proceed as follows; The potential of a point charge q 
located at 4^, and facii^ an earthed cone (ju^ = Constant) alone (i. e. , at first, 
without a spike and sphere) would be 


Vie = E - «’o> 


for 7 j < tJq, and 

oo 

Voe = E •"<■» - <Po> (AlO) 


for rj ^ tJq. It may be verified, by inspection of equations (A9) and (AlO), that V-^ = 
at ri = ri^, and that both solutions are finite everywhere and symmetric about the 
charge point | , ?? cp . Also, V. = V. . = 0 on the cone because of our choice of 

) = 0 . 

n '^c' 

Suppose now that the charge q is placed between an earthed sphere of radius R 
and the cone with an ellipsoidal spike attached at the cone’s apex (fig. (9)). The contour 
of the spike is given by 7 ] - rig = Constant. We must now superimpose on the potential 
of the cone another potential that will be zero on the cone and will give a resultant zero 
potential when v = r]g on the spike and t] = ri^ on the sphere R. For the spheroid 
with the semi -axes c and b. 





c 



and, from equations (A2) and (A3), 



+ 1-4 


2 



(All) 


(A12) 


to a high degree independent of 4, because 4<1 R/c»l. Notice that, since a 
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spheroidal coordinate system is used, the spheroidal surface V = Vg can be expressed 
by a sii^le value Tjg for the entire range of the other coordinates, | and <p. But two 
coordinates, 77 and | are needed (eq. (A 12 )) to describe exactly a spherical surface 
(R = Constant). Since now in figure 6, both 77 = 1 and 77 = «> are excluded, the super- 
imposed potential Vgj must have the linear combination Bj^j^P^(77)+ 

V^l = £ f A„„pS'(-{)[b„„C(^) +C^„Q|;'(,)]cos (A13) 

When equation (A 13 ) is added to equation (A 9 ) taken at 77 = TTg, the result is zero; 

<'^ic + = 0 = + c^„q;;*(^) (AH) 


And when 77 = 77^^ on the sphere, + Vgj is zero: 


(A 15 ) 


When equations (A 14 ) and (A 15 ) are solved for and 


Bmn = Qn <%) 


m... . C(^o)Qn - C(^S<(^o) 




^m 


m. 


.m. 


(A 16 ) 


Cxnn = -Pn 




K(^sK(nR) - K(vrY^(vs) 


(A 17 ) 


Adding this new potential Vgj to yields the complete potential Vj, when 77g < 77 < 
^o 


Vj = EE A, 

n m=0 


mn 


r(-«» { [ Vn * <'M C*") * 


And similarly, when Vo< V< Vr 
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OO 

v„ - E 5) h. [p”(,„) + Jq,”(„) I cos M<p - ?>o) (A19) 


To determine the constant a new variable <p' = cp - cp^ is introduced and an ex- 

pression for the flux density originating in the surface r] = of the charge is written. 
This expression is evaluated by invoking Gauss' theorem applied to a box fitting closely 
an element of the spheroid 77 ^. Using n for a "normal” (while n denotes the harmonic 
number), 


a_ = 
n 





(A20) 


Here, h^, h 2 , and hg are the metric coefficients (distance factors) appropriate to the 
orthogonal curvilinear coordinates of our spheroids. Specifically, 


h 


1 



hg = c 


2 



(A21a) 


(A21b) 


hg = r = C 2 


^(1 - - 1 ) 


(A21c) 


Returning now to equation (A20), the indicated differentiation may be completed with the 
help of equations (A21b), (A18), and (A19): 



2 

EE A 

n m =0 


mn' 


=“(-{) cos „)?;“(,) - Q;™(^)p“(,o)],^ 


(A22) 


where the prime on P and Q indicates differentiation d/dt]. Since our charge q at 
^o’ is a point charge, the function assumes in the limit the 

constant value ° 
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Multiplying the left side of equation (A22) by (-^) cos p<p’-hj^-hg x d| dcp' and in- 
tegrating over the surface of the spheroid within the cone at r] = gives 


/ 2it ^ fZir /*?c 

/ o_ ' P^(-|) cos p(p’h.h« d| d^’ = / / a (-f) cos p<»' • dS (A23) 

^ •= 0 *' 4 =-l n s 1 J ♦'O •'-1 n " 


• 2 jt r^c 


Now 



is continuous across the spheroid except over the infinitesimal area dS^ at cp = cp^ 
or <;<?' = 0 (cos p^' = 1), where the charge is located. Thus 9Vj/3n = 9V^/9n except 
over dS^, and the integral vanishes elsewhere. Therefore, by Gauss’ flux theorem 


2v 4c __ 

/ / P^(-^) cos p<p’ . a_ • dS = q • pf(-4^) 

0-1 ” 


(A24) 


Now multiplying the right side of equation (A22) by the same factor as before yields 



(A25) 


because the integral over dcp' is zero unless p h. m and the integral over d4 is also 
zero for s n. But a fundamental relation of spherical harmonics gives (ref. 6, 
ch. 5.212) 


I 
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(A26) 


<(voK^(vo) - qT(voK(%) = 


(- 1 )™ (n + m)r 
(? 7 q - l)(n - m)! 


where for noninteger values of n 


(n + m): = r e^* dt 

0 

> 

oo 

(n-m)!=/ t(^-”")e-tdt 
0 

and 

3| 3n 



Thus, from equations (A24) to (A26), and (A28) 


A 


mn 


(2 - 8^)(-ir 

2776 ^ 02(1 - 4) 



(n - m)'. 
(n + m)'. 



( 2 n + 1 ) 


(A27) 


(A28) 


(A29) 


We are now prepared to write the expression for Green's function G(P, P ), P and P 

O’ o 

having the coordinates n, <p and 77 ^, <p^ of the reference point and source point, 
respectively. Putting q = l (unit charge) into equation (A29), we get from equations 
(A18) and (A19) 


Po> = ZZ AmnC<-«{[Smn + ^ (^30) 

n m =0 

for 7?g ^ 77 ^ 7 ?q, and 
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Go<P.Po> =ZZ * [Cmn * P”(»o>Kw}“® M<t>-<P^) (A31) 


n m=0 


for Vo< Vr 


For a given space -charge density distribution PgCPjj) inside the sphere R and a 
given potential distribution V = Vg on the surface of the sphere, the exact expression 
for the potential V(P) at any point P located inside the space bounded by the sphere, 
cones, and spike is given by equation (Al); 


V(P) = / C 
«/vo! 


G(P, Pg) 
volume 




/ ^s<Ps>(^) 

t/surface 


dS 


(Al) 


In the volume integral the integration needs to be extended only over regions where 
Pg 7 ^ 0, that is, where trajectories exist. The second integral requires the knowledge of 
the normal derivative of at the spherical surface. In analogy to equation (A20), 





OO 

* [Cflo) * cos m(,^ - 


n m=0 


X(A32) 


P- 

d-q 




Q' 

dq 


r?=% 


Substituting equation (A29) for A^^^, and equations (A16) and (A17) for and 

respectively, and carrying out the indicated differentiations in equation (A32) yields 
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(2 - s;;^)(2n + 1) 


277€^(c2-b2)(l-4^)(4-4^) (%-l) 


1/2 


4/2 


^ ^ - C( ^)Qn (^q) 


n m=0 


( 


cos m{(p - cp^) 


(A33) 


ap“(-^) 8P“(-^)' 


9^ 9n 




Sometimes it is advantageous not to have the spike, if, for example, the potential 
depression and velocity -sorting due to a large space charge itself are stronger than de- 
sired. In this case, there is no need to use expansions in elliptical harmonics, but 
rather spherical functions are appropriate. We still use noninteger orders n of p“(-m) 
such that, on the cone, = 0* We forego the now much simpler derivation of the 

Green's function for a cone 9„ surrounded by a sphere of radius R and write the ex- 
pressionsfor G°, G°, and (9 Gq/9*)*=r (the superscript o denoting the absence of 
the spike). 


oo 



n m=0 




© 




(2n+l) 


) 


cos m((p - (Pn) 


(A34) 





= -(2n+l) 


n+1 

♦o 


r «+2 


E E 

n m=0 


cos m{cp - (Pq) 


(A36) 
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2 - 6 


m 


n 


mn 




3P“(-M) 3P“(-M) 


3fi 


3n 


M=M, 


(A37) 
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APPENDIX B 


TRAJECTORY EQUATIONS FROM LAGRANGIAN FORMULATION 

The basic geometrical configuration of the tube -collector arrangement is shown in 
figure 10. After leaving the last interaction gap, the spent electron stream passes 
through the fringes of the focusing magnetic field, if any, which leaks through the 
2a opening in an ideal (magnetic) shield into the collector region. Although sperical 
(or spheroidal) coordinates are most suitable for the trajectories inside, cylindrical 
coordinates lend themselves best for the transition region, which consists of a cylindri- 
cal tube terminated by the shielding plane which merges into the collector region. We 
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shall develop first the trajectory equations in spherical coordinates 1 , and <o; ') and c 
being the polar angle and the azimuth, respectively. The magnetic field is axisymnietric 
and derivable from the only components z) of the vector potential A, or a cor- 
responding scalar potential The expression for the generalized Lagrangian L* of 

a particle with a charge q is 


L* = T- q- V + qiT-A 


and for an electron with a charge q = -e 

L* = T + eV - eu ■ A 

Hereinafter e shall indicate the amount of the electronic charge e = +1.6x10”^^ 
coulombs, the sign having been taken care of. In our case of spherical coordinates, 

L* = -mfi^ + -*^ + i^(sin^ + eV - e *(p(sin ‘))A^,,(.^, 3) (Bl) 

2 

In equation (Bl), A^^(.,», 3) is assumed to have been transformed from A^^(r,z) by means 
of the relations shown in figure 10: 


r = 'jt sin 3 


(B2) 


Z ::: R - OB = R + '* cos 3 


(B3) 


dr^ = d*^ + -»^(d3)^ = d-*^ 



(B4) 


Here, r, cp, z are the cylindrical coordinates, and «, 3. <p the spherical coordinates. 
The momentum equations are obtained from 


A (v^L*) - V,l’ = 0 (B5) 

where q and q are the coordinates and the rate of change of coordinates with time, 
respectively. To evaluate equation (B5) we need to know the expression A = dA/dt, 
which appears in some components of equation (B5). Since A = A 3) has only one 
component A^ and ?A^ = 0, 
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dA „ 3A „ 
A =—^ = —£. 
^ dt 8t 


■I 


8A^ 

— ^ 

aqt 3* 


8A „ . 0A ^ 

_^+ 3 ^ 

3-9 


(B6) 


The application of the Lagiangian routine to equation (Bl) yields the following components 
for the equations of motion: 

3A 


or 


- -i(sin^ -'>9^ - 7 ] 


(sin^ 9)(^^ - v9^ - 7 j J — - ♦(sin 9)^ [■ “ ® 


~ - A (sin 9)<p - v(sin 3)cp 

? ^ 


3V 


<P 


dv 


= 0 


3 V 


L (B7a) 


2v*9 + ♦^9 - sin 9(cos 3)cp^ - 77 


— - ♦(cos '9)^A - r(sin 9 )^ — 


3A 


39 


<P 


39 


. . 2 • ’ 2 '2 
= 2'»'a 9 + * 9 - ♦ sin 9(cos 3)cp - rj^ 


— + T7p'#^(sin 9)^ • B,j(-*, 9) 
39 ® 


2 ^'i(sin 9)(J? +*(cos j + ♦(sin 9)^’ + 


~ + (cos 9)9 + — 

♦ A , 

(P. 


- 0 


= 0 (B7b) 

(B7c) 


In place of equation (B7c), we shall use Busch's theorem for determining Ip. From 
Busch's theorem for axisymmetric magnetic fields, 


27rr^ 


(B8) 


;//(r. z) = 277 



(B9) 


and i// (r , z ) represents the initial magnetic flux of the electron at r^, z^ at the time 
of its start from the cathode surface. In equation (B8) cylindrical coordinates were used 
for <p and the fluxes \p in agreement with our previous scheme. The expressions for 
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A and B = curl A = -grad in both coordinate systems are given in appendix D, 

ip lu 

together with the exact solution for an axisymmetric magnetic field bounded by an ideal 


shield with an opening hole . 

The use of equation (B 8 ) permits the elimination of <p from equations (B7). To 
eliminate the time dependence, the expression for the velocity u( 4 , 8 , ip) is used in the 
collector region, assuming, as before, steady -state case: 


u 


2 



dt 




(BIO) 


and 


^ ^ (Bii) 

dt 53 dip dn 


because d*/dip = 0. Substitution of equation (BIJ) into equation (BIO) yields 


u^(-t, 9, (p) = (♦’9)^ + (*9)^ + (4 sin 9 cp)^ = ->^9^ 


[■ •(«■] 


+ sin^ 9 


= 4292 




47t 2 4^ sin2 9 


1 [^(4.3)-i^J^ (B12) 


Now, for an electron entering the opening hole at a position r = R and 9 = 9. with a 

p _ / ■*' 

velocity Uj( 9 .), its velocity u (4, 9, 0 ) will be (writing 4 = pR, p = ♦/R): 


u2(p, 9, ip) = rV^^ 




vl [4'ip, 3) - 


.22 -2a 

4 t 7 4 sin 9 


= Uf(9.) + 2n^V{p, 9) 


(B13) 


2 

Equation (B13) may be normalized in terms of the dc electron energy Bivid- 

ing both sides of equation (B13) by u^ = 


2 

u2(p,3, cp) _ VV ^ V(P, 9) _ ^ ( 3 .) ^ V(P, 3) 

u2 u2 


(B14) 


or 
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(B15) 










V(P,3) . 


(47r2)(2r?gV^)(RV sin^ 0) 


Solving equation (B15) for 


yields 



vl [Kp, '5) - ^ 


47t^rV( 4 sin^ 3) 



(B16) 


It is possible now to eliminate the time from equation (B7a) such as to yield a time-free 

1 f 

trajectory p = p(3, 3., p p.(3^) and Pj^(3p being the initial values at the beginning 
of integration for a given initial energy class e^. The validity of the "time -free” tra- 
jectories has been discussed in the section Initial Beam Conditions. Here, we wish to 
repeat that this assumption holds exactly only at a given phase of entry (time) of the 
radiofrequency field into the collector and represents a good approximation at the mean 
phase of entry of the decelerated electron bunch. Since, in well-designed tubes, the elec- 
tron bunch accounts for about 80 percent of all electrons during a full radiofrequency cycle, 
our assumption describes a good average. With this understanding in mind we write 
= ■i'S; 'i' - d^/d-9 and substitute this into equations (B7a) and (B7b) to 

obtain 


p"-9 


2 



- p (sin 


2 




?A 

(sin 9)A ^ - ^ p(sin 9)<p — ^ - 0 
0p R R dp 


(B17) 


3 = -2— 3 + sin 3(cos 3)^^ + 


% ^ 
p V 


7p 


% 


+ — (cos 3)(pA^p + — ™ sin 9(p 


pR 




pR 


as 


(B18) 


After substituting 3 from equation (B18) and 3 from equation (B16) into equation (B17), 
the latter can be solved for p = p(3) providing that p-(3.) and pj(3j) are known for each 
single electron. 
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APPENDIX C 


DERIVATION OF CONTINUITY EQUATION 

We derive now an expression for the charge density in the collector. To make the 
problem amenable to numerical treatments, which we wish to limit to a reasonable effort, 
we invoke the following assumptions made earlier in this report (see section Initial Beam 
Conditions). 

(1) The time dependence is dropped; that is, all electrons entering the collector 
durii^ the interval of one radiofrequency cycle are thought of as entering at the same 
time (steady -state solution). 

(2) The current density across the beam is constant. 

(3) The axial velocity (energy) distribution is given quantitatively by large signal 
computer results or other accurate methods. 

(4) The angles of injection 0 !^^ are zero on the axis and proportional to the radius of 
injection away from the axis. 

(5) Compared to the size of the collector, the entrance area in the injection plane is 
very small (but not zero); that is, it can be regarded as approaching a point source. 

The last assumption is very important for achieving a numerical solution com- 
patible with present knowledge and an acceptable computer volume. It is that after a 
distance large compared to hole radius a but small compared to collector radius 
R » a, the electron trajectories are, in essence, independent of their initial radial 
position (point source) and, therefore, are determined only by the magnitude and direc- 
tion of the initial velocity vector. Thus, it is permissible to locate the slowest electrons 
at the edge of the hole and the faster ones more toward the center as long as their popu- 
lation and direction are kept in accord with the actual distribution function. 

To solve the continuity equation we can use either its differential, steady -state form 

V • (Pg • u) 0 (Cl) 

or its integral (valid generally also for nonstationary flows) 

Pg(xj) • J - , g(a.) i=l, 2, - (C2) 


where 


a(Xj^, X2, Xg) 

0(aj, ag, ag) 


(C3) 
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is the Jacobian of the spaces d(x) and d(a) occupied by the compressible electron fluid 
at the times t and t = 0, respectively. To be valid in the general case of motion of the 
electron fluid, including overtaking and crossovers (nonlaminarity), the integration 
should be performed "at large" in the Lagrangian fashion by tracking all electrons from 
their startup volume element d(a) at the time t = 0 to their respective position d(x) 
at t = t. We shall solve the problem for two different cases. In the first case we as- 
sume either that the beam has been refocused or that its radial velocities are so small, 
even without refocusing, that it may be considered essentially parallel, say of less than 
5° angular spread. 


Beam Almost Parallel 

Consider now figure 11(a) showing a tube of electron fluid injected at the hole be- 
tween and S„, at distances r. and r. + dr., respectively, away from the 
axis 0 - O' . Since the total dc current into the collector is continuous and the current 


n 




(a) Nearly parallel beam, (b) Diverging beam 

Figure 11. - Derivation of continuity equation. 
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density is constant across the beam, the current element dl^ crossir^ between 
and S2 is 


To compute the trajectories of the electrons which started between and at 
t = 0, we need to know their initial velocity vector Uj and also the injection radius as a 
function of energy. Suppose now that the axial energy distribution function is given by 
N(e): 


/ ^max 

= N(<^max) -N(«min) = = "1 (C6) 

^min 

because = 0. hi equation (C5), N(e) determines the fraction of electrons with 

energies larger than e. Since the distribution is non-Maxwellian, e is finite and, 
in most cases, — 2; is slightly larger than zero; and most of the electrons 

are grouped around e = 0. 3 to 0. 7. 

We return now to equation (C4). In light of what was said before regarding the dis- 
tribution of velocities, let the area dA(ej) in the injection plane occupied by electrons 
in the energy interval between be 

dA(ej) - 27rrj(e.) dr.(ej) (C6) 


On the other hand, the number of electrons per second (current) in this energy group is 
dl(€j^). This gives a current density 


dl(€.) 

= J_ = Constant 

27rr.(€j) dr.(e.) 


(C7) 


Equation (C7) permits the calculation of 


2i7 


♦=r. e. 

r 1 r ' N(^i) - 1 2 

I r(€j) dr(cj) = — I dl(ep = = 77 r?(e.) 

*'r =0 K ■ "^o 

^ mm 


After a flight time t, the infinitesimal volume element dA^^ x dh has moved along the 
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trajectory to CD in figure 11(a). Let the particle velocity at CD in the direction 
normal to CD be u , , then the height of the volume element at AD is dh • u /u. . 
From figure 11(a) it can be found that 


CD = 


P dp 

r -2 2I 

[p'^ + (p'r 


R 


(C8) 


Thus the initial volume element dAj^ x dh has changed into 


27tR 


2 p dp sin 9 


dh — 

■ 2 ^i 

p-2 + (p')2J " 


(C9) 


From equation (C2) we get, since the volume element as it moves along the trajectory 
contains an unchanging number of electrons: 


. dhp. = 27rR^pdpsin 8 ^ 

11,- .,1/0 .. 


dA 




1/2 




ri(6i) ?r. 


Pg(P 3) = (P^Uj) 


(CIO) 


. Vp^ + (p'f _ ^ P dp yjp^ + W'f (CII) 

Uj^(p,3) R^pdpsinS r\ (p,d) 


From equation (B14) we have 


u^(p, :9) = 




(B14) 


However, from equation (B 15) the velocity component u normal to CD is 


9 I V(o 01 ^ ^ *^o1 

^o (47r^)(277gV^)R^r sin^ 9 


(B15) 


because the rotation in the (^-direction does not contribute to the outflow of the electrons. 

It is obvious that the Jacobian (eq. (C3)) of our problem has been derived from geo- 
metrical considerations, rather than formally from the Jacobian determinant. As dis- 
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cussed earlier, for the solution (eq. (Cll)) to be valid generally for nonlaminar flows and 
non -steady -state conditions, the trajectory integration must be carried out for a suf- 
ficient number of initial times and all trajectories must be computed continuously. If 
crossovers and overtaking occur; that is, if streamline tubes which start from different 
initial segments of the injection plane or streamlines which take off from the same ini- 
tial segment but at a different time t^ are found at the same time t in the same volume 
element, all particles at crossover points must be added to account for proper space - 
charge density From this fact, the numerical complexity of integrating nonlaminar 

flow is apparent. With this remark we proceed now to discuss iionparallel, multivelocity 
streams. 


Diverging, Multivelocity Beam 

To treat this more complicated case, we need a two-dimensional distribution func- 
tion which relates energy and direction. From assumption 4 (at the beginning of this 
appendix), which is quite accurate for confined -flow beams emerging through an aperture 
in a magnetic shield, we have for the angle ^ Ib^) 


tan ^ 



(C12) 


The constant c is determined from 




(C13) 


C - 


“tmax 

r. 

1 , max 


Consider now a small element of the injection area (see fig. 11(b)) 


dAj = 2ttt^ dr^ 


(C14) 


The current injected through this area dl(rj^) contains electrons covering the entire 
energy spectrum - ^i - ^max s ame angular direction tan = cr^. 

Thus 
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(CIS) 


dI(r,)=J„. dAj 

The fraction of electrons per second entering between r. and r^ + dr. and with an 
energy between and + de^ is, therefore 

dl(r.) = M de. = J (2ffr.)^ (cos a) dr. de. (C16) 

de. ^ ° ^ dej ^ ^ 

The initial volume element containing a certain amount of fluid (electrons) of a given 
velocity class u^^j is 


2ttt^ dr^ cos Oj dh (C17) 

The charge density p.^. associated with the velocity class u^^ times the volume element 
(eq. (C17)) gives a charge 


27rr.(cos a.)p.- dh dr. 

1 '' ij 1 


(C18) 


As this charge moves along the trajectory, the volume element occupied by it will be 


2ttR^ P ^ 


dh 


} 


‘^ii 


(C19) 


and the charge density in it is p ..(p, 9 ). Thus, from conservation of charge we have 


Peij<P' " 


?p f \^3/ 


(sin 3)R^u.j^(p, 3) 


- — cos 0!i(PijU.p 


(C20) 


This is the space -charge density at the location p, 3 due to electrons injected at r = r . 
with an angle and within an energy class Cj. Because the motion is nonlaminar, there 
may be contributions to p ..{p, 3 ) stemming from electrons which have started at dif- 

eij 

ference initial locations r^ and with different energies Cj. For a steady-state motion 
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the crossovers will occur at certain time -independent locations. Thus, the total 
Pg(p, 3) is obtained, whenever crossovers occur, from the summation 



(C21) 


In equation (C20) the product cos represents the component J^j of the current 

density due to the velocity class j. Integration over all velocity classes results in the 
current density h J^, which is, by assumption, independent of location r. and, 
therefore, constant over the injection hole. 

Because trajectories are not known in advance, a closed-form solution is not 
feasible. Also, due to the interdependence of trajectories and space-charge densities, 
a relaxation -type solution is required. The reduction of space-charge density, which is 
desirable for reasons described earlier (in the section Space -Charge Effects), substan- 
tially speeds up the convergence. 
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APPENDIX D 


SOLUTION OF THE MAGNETIC -BOUNDARY-VALUE PROBLEM 


Expression for Component A^(r,z) of Magnetic Vector Potential 
in Cylindrical, Axisymmetric Coordinates (ref. 6) 


A^ (r,z) 


ZJ-1)^ ^^B(0,z) /r\(2n+^) 

n’.(n+l): U/ 

n-0 


= B(0,z) 


B 


1:2: 


o3 2’. 3; o5 


B 


3’, 4: 




(Dl) 


Expression for Leakage of Magnetic Induction Density Through a Circular 
Aperture in a Shielding, Infinite Plane Sheet 

Let the location of the plane in a r, c/?, z coordinate system be at z = 0, the plane 
extending to infinity in every direction normal to the z-axis. Let the sheet be made of a 
magnetic material with a ji - 00 and a B at saturation Let the B to the 

left of z = 0, that is, at negative values of z, approach a constant value B^ such that 
B(0, - -x;) = Bq. To the right of the shielding plane, B(0,z) decays to zero such that 
B(0, + 00 ) = 0. The scalar magnetic potential V^^ from which the required leakage field 
can be derived as B = -grad V^^ is (ref. 6, ch. 5.272); 

V^ = aBo^[(,-S) + i & cof\-f) + (D2) 

In equation (D2). ^ and J are the coordinates of oblate spheroids. The third coordinate 
is, of course, the angle w, same as in cylindrical coordinates. 

The coordinates | and t are related to circular cylindrical coordinates r and z 
through the relations 


z = a^f 


(D3a) 
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r = aj^l + 


)(1 - 




(D3b) 


The range of | and g is 0 < ^ < +1 and -«><£< +<=o; also, sign ^ = sign z. On the 
sheet, 1 = 0; inside the hole in the plane, t = 0; at the edge of the hole r = a. z = 0, 
both I and ^ are zero. Equations (D3) can be solved for ^ and 




2 2 2 -J 2 2.^ 2, 2 o 2 o 2^ 

: -a +z + T,r +z) +a(a -2r +2z) 


(D4) 




(D5) 


2 2 2 J, 2 2x2 2/2 o2o2 

r - a + z + Jr(r + z ) + a (a - 2r + 2z 


Notice that inside the hole in the plane z = 0, f 0. The values ? in the hole may be 
obtained from equation (D3b) by setting f = 0 and r = a(l - At r = z = 0, it 

follows 1 = 1 from either equation (D3b) or equation (D5) as the limit of | for 
z ^ 0 (r = 0). The transformation from cylindrical coordinates r, z into spherical 
coordinates ->, 9 may be accomplished with the use of equations (B2) and (B3); 


r = * sin 9 


(B2) 


z = R + ♦ cos 9 (B3) 

Substituting first equations (B2) and (B3) for r and z into equations (D4) and (D5) and 
then the values ^ .9) and | = ^(*, 9) into equation (D2), an expression 

9) is obtained, from which we get 


B,U,9) = - 


av. 


m 




1 av 

BsU,a) = -i “ 


^ dS 


(D6) 


(D7) 


This procedure is straightforward. Although the resulting equations can be obtained in 
closed analytical expressions, they are very cumbersome and are more easily produced 
through numerical analysis. 
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Equation (D2) and the field components (eqs. (D6) and (D7)) are valid everywhere. 
Inside the hole r < a, a series solution for A^(r,z) may be obtained from equations (Dl) 
and (D2). Now on the axis r = 0, § = 1 and z = af. Thus, equation (D2) simplifies to 


and 




(D8) 


B^(0,z) = 


3V^(0,z) 


= 


dz 


1-icot-l^^ 


1 B 


o a 


TT /„x2 

1 + 


f) 


(D9) 


Notice that B (0, 0) = — B , which is the same as the field of a long homogeneously wound 

Z Q O 

solenoid at its end. 


Since B = curl A, we get 


„ , . aA (r,z) 

B (r,z) = + 

r 0r 

3A^(r,z) 


Bj.(r,z) = - 


dz 


(DIO) 

(Dll) 


After B^(0, z) (eq. (D9)) has been substituted into equation (Dl), expressions for 
B 2 (i’, z) and B^(r, z) are obtained from equations (DIO) and (Dll). The first six deriva- 
tives of B (0, z) are listed below. Using j? as an abbreviation for z/a. 


/ 1 \ 2B 

B^^^O.z) = P 


(1 + 


(D12) 


B^^^O, z) = ^ 


77a^ (1 + j?2)3 


(D13) 


b(3)(0,z) = .!!^1:^- 

77a^ (1 +j?^)‘^ 


(D14) 


Bf(0,z).--^^. 

/i . ., 2 \^ 


48B, 


3 - 


(1 + ^^) 


(D15) 
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(D16) 


B<5)(0, z ) = * 35^^ 




2.6 




2 4 


^ ^ . 9. * 


iia 


(1 


From equations (D12) to (D17) and (Dl), 


A^^(r.z) 

B„ 


1 - i cot 




27ra2 (1 + ^2)3 


(D17) 


- 3-_5^. _ 5rJ_ ^ 3 - 14^2 + ,.4 ^ 

Swa"^ (1 + ^^2)6 48n-a^ (1 + 

For the equations of motion in appendix B this expression for A^^,(r,z) must be trans- 
formed into spherical coordinates P and 9 as before. Notice that z = 0. r = 0 (center 
of the aperture in fig. 10) does not coincide with -k = Rp = 0 (center of the sphere and 
origin of the spherical coordinate system). 

From equations (D18), (B2), and (B3), we get for = A^^,(p, 9): 



H 

(D19) 
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APPENDIX E 


COMPUTATION OF RESULTS 

The electron trajectories and equations associated with their derivations were com- 
puted on a high-speed digital computer. The details of the computations are not given 
here, but will be dealt with in a separate publication. Here we state the chosen approach 
only. 

The solution of the boundary -value problem was accomplished in part with the 

Green’s function series on the sphere and in part with the finite difference method along 

the spike since working with elliptical harmonics proved more involved than a direct 

way of forcing the boundary value along the spike. For the solution concerning conical 

boundaries, noninteger -order p Legendre polynomials P^(- cos >9) were needed which 

are zero along the cone (cos 9 = Constant) or hyperboloid (f ^ Constant) tangential to 

ni ^ 

the cone. Since the polynomials P^ (z) may be expressed by hypergeometric functions 
which contain the orders p and m as parameters, the latter may be computed from 
suitable hypergeometric series F (m - m + j/ + 1, m + 1, sin 9/2). See, for example, 
references 7 to 9, from which P^(cos 9) is obtained as 

Voices 9) = F 

2^(p - m)!m! 

Here the unknown is the order p which, for a given value of m and for cos 9 = Constant, 
makes P™(cos 9) = 0. 

A particularly difficult computation problem is posed by the requirement to solve the 
Poisson equation for current densities approaching 100 amperes per square centimeter 
at the injection hole. As many as six iterations were necessary for achieving a con- 
vergent solution. Since with every single iteration the results tend to "swing out" to 
far off, the application of the predictor -corrector method proved necessary and suffi- 
cient for speeding up and obtaining convergence. Variable-size mesh was used because 
of the congestion of trajectories in the vicinity of the injection hole. The computational 
step in computing the "high-space -charge" cases of trajectories was 0. OOIR and 0. 001° 
in 9 as minimum and 0. OIR and 0. 1° 9 as maximum. 




m -I- p + 1, m + 1, sin 


•■‘S 
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APPENDIX F 


SYMBOLS 


A(r,z) 

vector potential 

n 

"^nm’^nm’ ^nm 

expansion coefficients 


a 

radius of hole 

P 

B 

magnetic field density. 



Wb/m^ 

Co 

b, c 

minor and major axis of 



spheroids of revolution 

Co 

e 

electron charge, Asec 


G 

Green’s function 

q 

h 

height 

R 

hj, h2> hg 

metric coefficients 

r 

I 

current 

r, (p,z 

h 

normalized radiofre- 

■4, 8, cp 


quency current 

S 

j 

Jacobian determinant 

t 

Jo 

2 

current density, A/cm 

u 

j 

v/rr 

Y 

k 

coupling factor; summa- 



tion number 

V 

L* 

generalized Lagrangian 


m 

upper order of Legendre 



polynomials 



electron mass 

e 


Oi, 

1 

N 

normalized distribution 



function 

m 

n, u 

lower order of Legendre 

£ 


polynomials 



normal direction point- 
ing out of integration 
area 

power, W 

Legendre functions of 
first kind 

Legendre functions of 
second kind 

electric charge, A - sec 

radius of sphere 

distance from axis 

cylindrical coordinates 

spherical coordinates 

surface boundary 

time, sec 

velocity, m/sec 

volume enclosed by 
surface S 

voltage 

norm.alized radiofre- 
quency voltage 

normalized length, z/a 

injection angle 

Kronecker delta 

normalized energy 

dielectric constant, 

8. 86x10"^^ C^/N/m^ 
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H(t?), E(0, ^icp) 

factors in product so- 

ct 

circuit 


lution of Laplace equa- 

e 

electron, electronic 


tions in prolate sphe- 
roidal coordinates 

i 

initial injection plane 

3 

cone angle 

1 

lens 

M 

cos '9 

m 

magnetic 


oblate spheroidal coordi- 

0 

source or dc 


nates 

ov 

overall 

l.r?, CP 

prolate spheroidal coordi- 

RF 

radiofrequency 


nates 

r 

radial 

P 

normalized radius, r/R 

sat 

saturation 


volume charge density, 
3 

si 

superimposed 

ri 

Ve - e/m^ 

Asec/m 

efficiency 

electron charge to mass 
r3,tio 

t 

Superscripts: 

tube 

total derivative with 

a 

surface -charge density, 
Asec/m^ 

T 

respect to time, d/dt 
partial derivative with 

4' 

magnetic flux, Wb 


respect to space, 3/ 9 

Subscripts: 

B 

c 

Brillouin 

collector 
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